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1. Introduction 


Due to constantly increasing systems requirements and complexity, de¬ 
velopments in decentralized and compositional analyses and control have 
received a lot of research attention. In most approaches to decentralized 
control, sufficient conditions for the behavioral satisfaction of each compo¬ 
nent or subsystem are derived. A large-scale optimization problem is then 
formulated for the global objective with the local behavioral sufficient con¬ 
ditions as constraints. The resulting solution algorithm generally requires a 
solution for the local problems in the loop. As a result, the approach does 
not scale well for systems with non-uniform[^ local analyses. In order to over¬ 
come this limitation, a unihed framework for describing system behavioral 
properties at both the subsystems and global level is inevitable. The use of 
IQCs seems appealing since they provide a way of describing relationships 
between processes evolving in a complex dynamical system, in a form that is 
convenient for analysis. 

IQCs were originally introduced as a way to characterize troublesome 
nonlinearities, time-variance, and uncertainties that plagued the application 
of linear control methods. Ever since, they have been used extensively in 
control PQ El E], and presently gaining interests in optimization la i , as 
a unihed way of characterizing system behavior from input-output point-of- 
view. They are used extensively in this work to describe desired input-output 
system characteristics. 

The following issues are addressed in this paper: 

Unified description of control synthesis problems 

This paper describes a way of looking at control synthesis problems as 
convex optimization problems involving parametrized integral quadratic con¬ 
straints. This exposes a paradigm that can help transfer results between 
optimization and control. Certain properties of the IQCs that are a certih- 
cate for the BIBO stability of the resulting closed loop system are provided. 
Moreover, the feasibility and the solution of the synthesis problem for linear 
systems are described. 


^this will make indexing very challenging - as the subsystems have to be handled 
individually and differently 
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Global admissibility condition 

Given a global synthesis problem associated with an interconnected sys¬ 
tem and a hnite collection of local synthesis problems associated with the 
local subsystems, the condition for which the feasibility of the local synthe¬ 
ses problems implies that of the global problem is derived. This condition is 
similar to existing results in literature ilEli, different only by the contain¬ 
ment of an optimized parameter. Solution approach using the alternating 
direction method of multipliers is given. Also, corollaries under the assump¬ 
tions of well-posedness and quadratic parametrization are also provided. 

Localization 

The term localization is introduced to describe the local syntheses prob¬ 
lems for which the global admissibility condition holds for all values of the 
optimized parameter. Furthermore, the concept of closest localization is in¬ 
troduced by including certain norm-like objective. It is then shown that the 
closest localization corresponds to the most relaxed local problem given any 
interconnected system and associated global problem. In addition, group 
localization is introduced as a way of localizing the global problem among 
disjointed sets of subsystems. An alternating minimization algorithm is also 
given, without convergence analysis/proof, as one way of solving the group 
localization problem. This relies on well established results regarding the 
convergence of alternating minimization algorithms [SlEniE]. 


2. Notation 


Throughout the paper, the following notations are used: M and M_|_ de¬ 
notes the set of real numbers and positive real numbers respectively. M” 
denotes the hrst orthant of the n-dimensional euclidean space. The set of all 
symmetric positive dehnite and positive semi-dehnite matrices are denoted 
by S++ and §+ respectively. The Euclidean norm of a vector x G M"’ is de¬ 
noted by ||x|| = (x'''x)^'^^. The quadratic form ||x||p = x^Px is dehned 
for any matrix P G S+. The expression P ^ Q means that the matrix 
Q — P G S+. The set of the eigenvalues of a square matrix A is denoted 
by A (A) with Xmax{A) = max{A(A)} and Amm(A) = min{A(A)}. Simi¬ 
larly, The set of the singular values of a matrix A is denoted by o'(A) with 
o'maxiA) = max{cr(A)} and crmm(A) = min{cr(A)}. Occasionally, the block 

for compactness purpose. 


symmetric matrix 


A 

B' 


C 


is written as 


W 

(^) 

A 

B 

* 

C 
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Let U r < n he a. full rank matrix, t/-*- denotes the orthogonal 

complement of U, i.e. Ulf-^ = 0 and [U f/-*-] is of maximal rank. 

The euclidean balls Bj.( 0 ) and Br(xo) are dehned respectively for some 
r G M+ as B^( 0 ) = {x : ||x|| < r} and B^(xo) = {x ; ||x — Xo|| < r}. 

The space of all square-integrable signals f; M+ —)■ M"' satisfying 

poo 

/ ||f(f)|pdf < oo 

Jo 

is denoted by L,^. Consequently, the £2-norm of a signal f G L2 is dehned as 

n2 = [r wmrdtY 

Let f : M+ h->• M"". Then for each T G M+, the function fp : M+ 1—)■ M" is 
dehned by 


Mt) = 


f(t), 0 < t < T 
0, t > T 


and is called the truncation of f to the interval [0, T], Consequently, the 
set L2e of all measurable signal f : M+ 1—)■ M"’ such that frit) £ IL2 for all 
T G [0, 00 ) is called the extension of L2 or the extended L2-space. The 
Fourier transform of f G L2 is denoted by 

poo 

i{juj) = / e~^‘^*i{t)dt. 

Jo 


3. Preliminaries 

This section develops some preliminary results required for the main de¬ 
velopments and applications in the subsequent sections. 

Definition 1 (Causality). Let Pt be a past projection operatoij^ dehned for 
any T > 0 by 

Prfit) = frit). ( 1 ) 

The operator A : L^e e-)- L^g is then said to be causal if P^A = P^APt for 
any T > 0 . 


^Operator is used to describe a function from one L 2 e space to another 
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Definition 2 (BIBCj^ Stability). A causal operator A : L^g e->■ L^g 
be stable if there exists c G M+ such that 

rT rT 


A{w(t)) A{w{t))dt < c / w{t) w{t)dt 


is said to 


( 2 ) 


for all T > 0 and w G L^g. 


3.1. Integral Quadratic Constraints (IQC) 

In general, IQCs give useful characterization of the structure of operators 
- which in this case are systems described by ordinary differential equations. 


Definition 3 (Integral Quadratic Constraints [2]). A bounded operator A : 
I^2e ^2e is said to satisfy the IQC defined by 11 if 



n(ja;) 


w{ju) 


doo > 0 


( 3 ) 


holds for all z = A(w), w G L 2 . 


Often, n is referred to as the multiplier that defines IQC and the short¬ 
hand notation A G IQC(n) is used to meean that A satisfies the IQC defined 
by n. 

Remark 1. This input-output characteristic will be used later to generalize 
the class of control synthesis problems considered in this paper. 

Remark 2. Observer that A G IQC(n) implies that A G IQC(an) for all 
a G M+. Thus, the set of multipliers /C(A) = {n| A G IQC(n)} forms a 
convex cone. This is an indication that control synthesis problems for a 
dynamical system described by the operator A can be described as corre¬ 
sponding optimization problems over the cone /C(A). 


Proposition 1 (Stability Certificate). The operator A : L^g >-->• L^g is stable 
if and only if there exist a boundec|^n(ja;) given by the conformal block 


n(ja;) 


nii(ja;) ni2(ja;) 

ni2(ja;)* 1122 (jw) 


( 4 ) 


satisfying nii(ja;) ^ 0, and Yl 22 {,i^) A 0 for all a; G M, such that 


A G IQC(n). 


( 5 ) 


^bounded-input bounded-ouput 
"^bounded eigenvalues 
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Proof. (^) Suppose A is stable. Then, there exists c G M+ such that 


/ A(w(t))'''A(w(t)) — cw{t)^w{t)dt < 0, 

Jo 

which, after using Perseval’s identity, implies that 

AelQc(^ 0^_°/ 

(«^) Let TTij A sup^ ||njj(ja ;)||2 , i,j G {1,2}. Then, for all to G M+, 


^w(ja;) 

A(w(ja;)) 


n(jca) 


^w(ja;) 

A(w(ja;)) 


du > 0 


implies, using the young’s inequality, that 




(VTII + e) ||w ||2 - 7122-— ||A(w )||2 > 0, 




for some e > 0. Choose e > —, and since w G L^g, it follows that 

7r22 


A(w(t))'''A(w(t))(it < c / w{t)'w{t)dt 


T, 


for all T > 0, where 


o<c4hii±hi 

71228 - 71(2 


( 6 ) 

(7) 

□ 


3.2. Dissipativity 

Consider a continuous-time, time-invariant dynamical system described 
by 


y . / x(t) = f(x(t),w(t)), f(0,0) = 0 

■ \ z(t) = h(x(t),w(t)), h(0,0) = 0 

with X G X C M"-, w G W C and z G Z C 
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Definition 4 (Supply Rate). Given the dynamical system in ([^, a supply 
rate is any mapping s : W x Z i—)■ M satisfying 



s(w(t), z(t))\dt < oo 


for all to,ti G M+ and for all input-output pair (w(t),z(t)) satisfying 


(9) 


(§ 


Henceforth explicit time arguments in signals will be dropped, except 
otherwise needed for clarity. 


Definition 5 (Dissipativity). The system E is said to be dissipative with 
respect to a supply rate s if there exists a differentiable and nonnegative 
function R : X i—)■ M+ such that 


VR(x)'''f(x, w) — s(w, z) < 0 


( 10 ) 


for all X G X, w G W. 


Remark 3. The inequality in (10) is referred to as the Dissipation Inequality 
and describes a property of the system that stipulates that within any time 
interval, the change in the internal stored energy cannot exceed the total 
externally supplied energy. Hence, there can be no internal “creation of 
energy”; only internal dissipation is possible. If (10) holds with equality for 
all X G X, w G W, then S is said to be lossless with respect to s. 


Remark 4. It can be seen that IQC generalize the dissipativity framework to 
supply rates that are themselves dynamical systems by allowing frequency- 
dependent requirements to be described. Consequently, a definition of IQC 
in terms of dissipativity can be given. In the line of the definition given in |6] , 
let {A, B, C, D) be a realization of a stable LTI system T with state vector 
7) G such that the multiplier H admits the factorization H = . 

Then the system described by (|^ is said to satisfy the IQC defined by H if 
there exists a mapping R : M"" x 


V,rR (x, 77)^f(x, w) -h Vry (x, r)Y 


-|- B 


{0x1 +D 


w 

z 



< 

Cr) + D 



( 11 ) 


for all X G X, w G W and z = h(x, w). Traditional dissipativity definition is 
recovered in the special case where 4/ is static. 
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4. Synthesis Problem 


Throughout the remainder of this paper, we shall limit our attention to 
static supply rates. Generalizing to dynamic supply rates is left for future 
work. If the static supply rate is parameter dependent, a fairly general control 
synthesis problem can then be dehned. 


u 

V 



Figure 1 : Dynamic Operator 


Consider the causal operator H : x e->■ as shown in Fig. 

Here, x G M""” represents the internal states of the operator, (v, y) G x 
L2g are the input-output pair over which the desired performance objectives 
are dehned and u : L^g x L^g e->• L^g is the control input which is to be 
designed to achieve the performance objectives. The numbers + Uy = 

n, Uu G N are the dimensions of the respective spaces. Now, we have all the 
ingredients to give a formal dehnition of a synthesis problem. 

Definition 6 (Synthesis problem). Given a parametrized static multiplieij^ 
X(7) G M”^"',Xii(7) ^ 0,X22(7) ^ 0 , with parameter 7 G M. The multiplier 
X (7) is monotonically increasing in the parameter 7, namely X (71)—X(72) ^ 
0 77 71—72 < 0 . The synthesis problem for the operator H : L^g xL^” L^g 
is given by the optimization problem 


J min 7^ 

\ s.t i/G IQC(X(7)), 


( 12 ) 


and is denoted by Ph(-^(7);7^) 

Remark 5 (Feasibility). The synthesis problem Ph'(X( 7);7^) is said to be 
feasible if the associated optimization problem is feasible, i.e there exists a 
control law u : x 1—)■ Ln“ and at least one parameter 7/- G M such 

that H G IQC(X(7j)). 


^the terms multiplier and supply rate are used interchangeably 










Remark 6 (Solution). The shorthand 7* = Pff(X(7);7^) is used to mean 
that 7* is the solution of the synthesis problem P//(X(7);7^) and is dehned 
as the solution of the corresponding optimization problem Ph in (12). 


Remark 7 (Properties). The following properties of the synthesis problem 
are obvious; 


• Pff(X(7);72) = PH{aX{'y)-,'y^) for all a G M+ 
. FH(X(7);7') = ^r/(X(7)+^^(7);7') 


Hence, without loss of generality, the supply rates (multipliers) will be as¬ 
sumed symmetric. The above properties also demonstrates that scaling is 
inconsequential for the solution of the synthesis problem. Only the eigen- 
structure of the multipliers are important. 


4.1. Basic Synthesis Problems 

Here, we discuss the synthesis problems corresponding to some basic con¬ 
trol objectives which can be described within the dissipativity framework. 
Consider the dynamical system 


f X = f(x,u,v), 
\y = h(x,u,v). 


( 13 ) 


The corresponding synthesis problem is given by 




mm 7 
s.t 


H ■ 


VH(x)Tf(x,u,v) 


V 

y 

T 

V(7) 

1 1 

'C <1 

1_1 

” and y = h 

(x, u 


< 0 , 


(14) 


4 - 1 . 1 . Square Systems 

The term square is used to describe systems with equally dimensioned 
input and output. First, the input-output objectives considered are dehned. 
Then, the synthesis problem encompassing the dehned objectives is given. 


Definition 7 (Positive-real or Passivity Objectives). The state space system 
P in ( 13 ) is: 


passive if it is dissipative with respect to the supply rate s(w, z) = w" z. 
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• strictly output passive if there exists e > 0 such that H is dissipative 
with respect to the supply rate s(w, z) = w'''z — e ||z||^. 

Definition 8 (Synthesis Problem for Positive-real Objectives). The synthesis 
problem for positive-real objectives is given by the optimization problem 


P, 


mm — £ 


s.t 


H ■ 


H G IQC 
£ > 0 , 


0 / 

/ -el 


(15) 


Remark 8. It is also possible to be interested in just the feasible problem. In 
which case, the objective of the synthesis problem is to certify a strict output 
passivity for a given e: > 0. The corresponding synthesis problem in this case 


is given by Ph 


0 / 

I -el 


;0 


4-1.2. Rectangular Systems 

Next, we consider systems in which input and output dimensions are not 
necessarily the same. 


Definition 9 (^ 2 -gain). The £ 2 -gain of the state space system H in (13), 
from the exogenous input v to the regulated output y, is given by 


7 (P) = inf {j\3b 3 Wy^Wl < 7^I|vt||2 + b, V v e VT > O} (16) 


Proposition 2 (Section 3.1 [ 12 ] )• The state space system S in (|^ has 
£ 2 -gain< 7 if it is dissipative with respect to the supply rate s(w, z) = 


Sometimes, it is desired to keep the system trajectory within a known set. 
That is, to enforce some constraints on the reachable set of the system. These 
constraints could translate directly to requirements such as safety, behavioral 
or validity. The following theorem describes the condition which provides a 
guarantee for the invariance of some sets under norm-boundedness assump¬ 
tions on the exogenous input. This type of objective is termed reachability 
because of the direct constraints imposed on the set of reachable states. 
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Proposition 3 (Reachability). Suppose f : X X W H- >■ M” in (|^ is Lipschitz 
continuous in both arguments. Given /3 > 0 and a differentiable R : X i—>■ M 
satisfying R(0) < (3. Suppose that ||w ||2 < {3. If 

VR(x)^f(x, w) < w^w (17) 


for all X G X and w G W, then V (x) <2(3 for all t G M+ - meaning that the 
set £ = {x| R(x) < 2(3} is invariant 


Proof. Integrating the dissipation inequality in (17) yields 


l^(x(f)) < R(0) + f ||w(r)|pdr 

Jo 

< 2(3. 


□ 

Definition 10 (Synthesis Problem for System-Norm Objectives). The syn¬ 
thesis problem for norm objectives is given by 

" min 7 ^ -f- \(3‘^ 
s.t 

p^. VR(x)Tf(x,u,v)- 

||yf-/3Mx)<0, 

^ for all X G X, V G and y = h(x, u, v), 

where A > 0 is a weighting on the multi-objective performance index. 

Remark 9. The first inequality constraint, together with the first term in 
the performance index, corresponds to the £ 2 -gain objective. The second 
inequality and the second term of the performance index define an H 2 -like 
objective for the synthesis problem. This is obvious when the system H is 
linear time invariant with a Lyapunov function V (x) = x^Px. 

Remark 10. Subsets of the above synthesis problem can be considered. For 
instance, removing the second inequality together with second term of the 

performance index results in the problem Ph 



7 ^/ 0 

0 -I 


< 0 , 


(18) 
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4-2. Linear Systems 

The goal is here to give a full description of the synthesis problem for 
linear systems. This includes the derivation of the feasibility conditions as 
well. Consider the linear system 



where x G is the state vector, y G M”*' is the output of interest, G M”"* 
is the measured output, v G M”” is the exogenous disturbance, and u G M”" 
is the control input. 

The controller considered is a hnite dimensional LTI system described as 



where Be, Cc, Dc are appropriately dimensioned parameters to be found. 
The controller input y^ and the output u matches the measured output and 
control input respectively of the LTI system H^. 

Consequently the closed loop system is given by 



Given a static multiplier X = X~^, it is straightforward to see that the 
condition Hd G IQC(X) is equivalent to the existence of a matrix P G S+ 
such that 
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The corresponding synthesis problem is then given, by including a parame¬ 
ter dependence on the multiplier, as The following remarks 

comment on the feasibility and the solution approaches for (X( 7 ); 7 ^). 

Remark 11 (Feasibility). Let 




Xn X12 

X 22 


It is well known ( see [13] Chap. 4, [T3| and references therein), using elim¬ 
ination lemma that Hd € IQC(X) if and only if there exists Qi,Q 2 £ 
such that the following hold 


Qi I 
I Q 2 


^ 0 , 


u 


T 



/ 



1 



Cl Pii 



where 


-1 T 


H 


/ 

' A Bi' 

T 



Cl Pii 



T 


0 0 

Ql 

i 

0 


1 

0 

0 

-X12 


Ql 0 

0 

0 


0 -XT 

0 

-X22 . 

- 

0 0 

Q2 

0 ■ 


1 

0 

0 

-X12 


Q2 0 

0 

0 


1 

0 

0 

-X22 J 



/ 

A Pi 
C-i Bn 


I 

A Pi 

C-i Pii 


B^O, 


(24) 


V^O, 


u = 


C2 

D21 


and V = 


B 2 

Di2 


Remark 12 (Solution). There are many ways to compute the solution of 
Ph-^;(X( 7 ); 7 ^). One simple method proceeds as follows; First, use bisection 
[T5| to minimize the sub-optimality level 7 . Suppose the minimum value 
of 7 is achieved by Ql and Q^. Compute non-singular Pi and P 2 with 
PiPj = / — Q 2 Q 1 , and determine P by solving 


Ql P2 
I 0 


P 


I 0 
Q2 Pi 


(25) 
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Once P has been determined, the condition Hd G IQC(X) is linear in the 
parameters Ac, -Be, Oc, Dc and a feasible point can easily be determined us¬ 
ing available efficient convex optimization routines. Other regularization 
constraints can also be added to improve the conditioning of the resulting 
controller, or even determine the order of the controller. See [13] for more 
details on the different ways to this computation. 

5. Problem Definition 



u, 


z ◄- 







M 


V.' -^ 




^ W 




-► y* 


(a) Global System (b) Local System 

Figure 2: Large-scale Interconnected System. 


Figure 2a shows the interconnected system under consideration. The 
schematic for each block is shown in Figure The interconnection 
output vector v constitutes the exogenous inputs Vj, i = 1,..., N for the 
subsystems. On the other hand, the controlled outputs y^, i = 1,..., N of 
all the subsystems constitute the interconnection input vector y. The static 

characterizes the interconnection 


matrix M = 
topology. The term 


'Mn 

Mi2 


y 


V 

_M21 

M22 


w 

z 


^‘global” 


is used to describe the entire interconnected 
system as shown in Fig[^ while “local” refers to the individual subsystem as 
shown in Fig[^ The goal of this section is to describe the synthesis problem 
for the global system in terms of the interconnection input-output pair (w, z) 
and the local systems in terms of the input-output pairs (vj, yt), i = 1,... ,N. 
Moreover, to establish a relationship between the local and global syntheses 
problems such that the solution of the local synthesis problem provides at 
least a feasibility certiheate for the global synthesis problem. 

Definition 11 (Global Problem). Consider the interconnected system G 
in Fig dehned by the local subsystems H^, i = 1,... ,N and the inter- 

{n^+n^)x(ny+n^) 


connection matrix M G 


Given a global multiplier iF(7), 


14 
























parametrized by 7 G M, the global problem is defined as the synthesis prob¬ 
lem Pg'(W^( 7 ); 7 ^) which is given by the optimization problem 


P, 


mm 7 
s.t 


G ■ 


w(jw) 


hh(7) 


w(jw) 

z(jw) 


(26) 


do; > 0, V w G L 2 ™ 


Remark 13. It is straightforward to see how the definition above can be 


extended to encompass the syntheses problems given in (15) and (18). The 
form above is used for simplicity of exposition. 


Definition 12 (Local Problem). For the interconnected system G, we define 
a local problem as finite collection of syntheses problems associated with the 
subsystems Hi, i = 1,N. In other words, given a finite collection of 

parametrized multipliers {Xj( 7 )}^]^, the local problem Pl (^{-^i( 7 )}^i ; 7 ^j 
is given by 

Pl ({V(7)}ii;7") = {Ph, (27) 


where each (Xj( 7 ); 7 ^) is the synthesis problem associated with each sub¬ 
system Hi. 

Remark 14. The local problem Pl ^{-^i( 7 )}^i ; 7 ^^ is feasible is if there ex¬ 
ists at least one 7 ^ G M such that all the syntheses problems P^-. (Xj( 7 j); 7 ^), 
i = l,...,iV are feasible. Moreover, 7 * = Pl ^{Xj( 7 )}^^; 7 ^^ is called 

the solution of the local problem Pl ^{-^i( 7 )}^i ; 7 ^j and is given by 7 * = 

max{ 7 j}^j^, where % = P^. (Xj( 7 ); 7 ^) is the solution of synthesis problem 
associated with the subsystem Hi. 

Definition 13 (Global Admissibility). The local problem P^ ^{Xj( 7 )}^^ ; 7 ^j 

is said to be globally admissible by the global problem Pg{W ( 7 ); 7 ^) (or sim¬ 
ply admissible) if the following condition holds for all parameter 7 G M; 



for all , w, and z. 



w 

z 


< 0 


( 28 ) 


15 
















Remark 15. If a local problem is admissible by the global problem Pq^ 
then the feasibility of Pl is a certihcate for the feasibility of Pq- 

Theorem 1 (Global Admissibility Condition). Given the static interconnec¬ 
tion matrix M G the local problem Pl ; 7 ^ 

with each Xj( 7 ) given by the conformal block, 


^*( 7 ) = 


xFh) xpi^) 
[xF'h) xph) 


(29) 


is admissible by the global problem P( 5 (hh( 7 ); 7 ^) if and only if the matrix 
inequality 


r ^“(7) 1x12(7) 


r Ml 

T 

— 1 X 22 ( 7 ) 



[Ml 

I 


( 7 ) 

X22(7) 


I 


-1^12(7) I -fGii(7). 


(30) 


where 


A'J*(7) = 


is feasible for all 7 G M. 


Af(7) 


xSb) 


, j,ke { 1 , 2 }, 


(31) 


Proof. The global admissibility condition in (28) is equivalent to the inequal¬ 
ity 


V 

T 

-A“(7) A'>"(7)- 

V 


w 

T r 

.y. 


A>277) A22(7)_ 

.y. 


z 



H"ii(7) ^^ 12 ( 7 ) 

»7(7) 1 ^ 22 ( 7 ) 



< 0 , 


holding true for all v = [v^ • • • y = [y 7 • • • y^]"'", w and z satisfying the 

interconnection constraint. After rearranging terms, it is straight forward to 
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see that the global admissibility condition is equivalent to 


V 

z 

T 


-iU22(7) 

X12(^) 

- 1172 ( 7 ) 


V 

z 

y 

w 


(7) 

-hFi2(7) 

X22(y) 

-hFn(7). 


y 

w 


V 

z 


= M 


y 

w 


for all V, z, y and w. 




r ^^^7) 

X^2(7) 

M' 

T 

— 11 / 22 ( 7 ) 

-W^S(7) 

I 


( 7 ) 

X22(7) 



L -iUi2(7) 

-VFii(7) 


□ 


Using the global admissibility condition, a relaxed optimization problem 
for the global synthesis problem can be written in terms of the local problem 
as follows 

Definition 14 (Relaxed Global Problem). Consider the interconnected sys¬ 
tem G in Fig dehned by the local subsystems Hi, i = and 

the interconnection matrix M G Given a global multiplier 

hF( 7 ), parametrized by 7 G M, the relaxed global problem is dehned as the 
synthesis problem which is given by the optimization problem 


/" . 2 

mm 7 

s.t 


Pg- 


M 

I 


T 


'11 


-iU22(7) 


' 12 ^ 


'12 


wUi) 


722 


-iUi2(7) 

^ 0,^22 ^ 0 , 


-iUn(7) 


M 

I 


^ 0 , 


(32) 
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Remark 16 (Solution). Similar to the approach used in |6], a solution method 
based on the alternating direction method of multipliers (ADMM) [16] is con¬ 
sidered. ADMM belongs to a category of operator splitting techniques. It 
solves convex optimization problems by breaking them into smaller pieces, 
each of which are then easier to handle. At each iterative step {k-th step), 
the solution implements the following Gauss-Seidel sweep: 


1. = argminPz where 

' N 

min7^ ++ 




r 

Zl2 1 

Ml 

T 

— 1 ^ 22 ( 7 ) 

-wUi) 

I 


Z 12 T 

Z22 



L - 1 x 12 ( 7 ) 

- 1 X 11 ( 7 ) J 


^ 0 . 


2. = argminPx where 


N 


mm 




in \\x, - zy + n 


^ IIF 


i=l 


S.t 


P, gIQC(Z,), z = 1,...,N. 


3. Dual Ascent Update 


(34) 


yk+\ 
* I 




k-\-l 


zy* 


(35) 


It is clear that step 2 is a candidate for parallel implementation on a dis¬ 
tributed processor. As a result, the ADMM solution steps described above 
can be interpreted as a method for decentralized control design. 
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Ultimately, we would like to translate the synthesis problem associated 
with the global problem to a corresponding set of syntheses problems asso¬ 
ciated with a globally admissible local problem. However, it is not clear how 
this can be done at this point. 

Moreover, the relaxed global problem in (32) incorporates the local syn¬ 
thesis problem. In order words, the solution of the global problem requires 
at least the feasibility conditions for the local problem. This makes it unsuit¬ 
able for cases where one is interested in designing performance specihcations 
for the local subsystems using only the interconnection and global objectives, 
without necessarily solving the local problem. 

These shortcomings are due to the general nature of the parameter de¬ 
pendence in the global admissibility condition above. It turns out that, by 
imposing a particular parametric structure on the global and local multipli¬ 
ers, the global admissibility condition can be made parameter-free. Hence, 
a much nicer decoupling of the global and local problems can be obtained. 
The following corollary examines the parametrization that is exploited in this 
paper. 


Corollary 1 (Well-posed Interconnection). Suppose the interconnection de- 


hned by M = 


Mil Mi2 
M21 M22 


y 

w 


I-)- 


are invertible, then the global admissibili 


is such that M 21 M 21 


Mf2X^\'j)Mi2-Wii{'y) 

My2X^\'y)-Wi2h)M2i' 

{Mf2X^\^)-Wi2{'j)M2iy 

X^\'j)-MlW22{l)M2l_ 


y condition in (|30|) is equivalent to 

^ 0 (36) 


Proof. Observe that since M 21 M 21 are invertible, the intercon¬ 

nection matrix can be factorized as 


Mil Mi2 


0 

Mi2 


/ M/]^M22 

M 21 M 22 


_M21 

0 


_M/2Mn I 


where M /2 


, -1 


= {MJ 2 M 12 ) MJ 2 and Mji = {M 2 iMj^y 
Pre- and post-multiplying the right hand side of the global admissibility 
condition in (30) by and P respectively, where 


P = 


I M|j^M22 
MI 2 M 11 I 

-1 




I 
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yields 



which after expanding and reordering the matrix elements becomes 


_ (MTXi 2 (^) - Wu{l)M2iy \x^\l) - MJ^W22 {i)M2i\ - 

□ 

Remark 17. The reqnirements that be invertible dehnes 

the well-posedness of the interconnection namely; the pair (v, y) is uniquely 
dehned by (w, z) and vice versa. 

Remark 18 (Local Passivability). Consider a passivity global objective, i.e 

W = ^ Q . It straight forward to see that, if the interconnection is well- 

posed as dehned above, the interconnected system is locally passivable if and 
only if there exists a block diagonal matrix D such that 

- M 21 = 0. (37) 

Equivalently, the interconnected system is locally passivable if and only if the 
matrix M 12 (M(^Mi 2 ) M 21 is block diagonal. 

Corollary 2 (Quadratic Parametrization). Suppose the global supply rate 
W ( 7 ) is quadratically depends on the parameter 7 as follows 

PP( 7 ) = + 27 IP 2 + IE 3 . (38) 

If the local supply rates are equivalently parametrized as 

^*( 7 ) = + 27 X ^2 -|- Xjg, (39) 

then given the static interconnection matrix M G the local 

problem (^{-^*( 7 )}^i j 7 ^) is admissible by the global problem PciW ( 7 ); 7 ^) 
if and only if the linear matrix inequality 

Q'IYlQi-QIYgQ2 P^, (40) 
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is feasible, where 


Qi = L 


2x2 


Mi 2 Mil 

0 I 


, Q 2 — I‘. 


2x2 


I 0 
M22 M21 


(41) 


Yl{X) = 




Xi, 



Xn, 


Xn 2 



X13 



Xn. 


Xns 


, Yg{W) = 


Wi 


IT/ 

ITs. 


(42) 


Proof. First, observe that the global admissibility condition is equivalent to 

^ 12 / 


1 _ 

T r 

0 / 



X22(7)_ 

/ 0 

M22 M21 


Mi2 Mil 

0 I 


hF(7) 


I 0 
M 22 M 21 


^ 0 , V7 e 


Next, substituting the parametrization in (38) and (39) yields 

T 


7/ 

I 


{QJYlQi-Q^YgQ2) 


I 


^0, V7 e 


^Q1YlQi-QIYgQ2Y^Q. 


□ 


Remark 19. Given the interconnection matrix M and the global supply 
rate parametrized as 

^ o,x“ ^ 

I N 


in (38), solving (40) for any feasible Yl satisfying 
0 ,X 22 ^ 


0 yields a globally admissible local problem 
Pl (^{^i(7)}ili ;7^j- This process can be interpreted as a supply rate allo¬ 
cation technique, in which performance specihcations are prescribed for each 
subsystem such that the global performance is guaranteed for any quasi- 
convex performance level. This structure is exploited in the next section to 
analyze the achievable performance levels by the set of all globally admissible 
local problems. 
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Remark 20. If the interconnection is well-posed as defined above, a more 
structured quadratic parametrization can be defined, viz. 


and 


W{-,) = 


V.(7) = 


7^11 + ITn 

27IT12 -|- IT12 

_( 27 lTi 2 + i^i2)^ 

1^22 

r 


_[2-jXy + xfY 

-1 


(43) 


(44) 


Thus the global admissible condition becomes 


M72XIIM12 - ITii 

0 

(M72X12 - W12M21 

- W12M21 

V 

* 

MJ^X^^Mi 2 - hhii 

M72X^' - IT12M21 

- IT12M21] 

T 

99 

X -MIW22M21 


^ 0 

(45) 


6. Localization 


The term localization is used to describe local problems obtained from a 
global problem in a way that certain global properties are ascertained. 


Definition 15 (Localization). Given an interconnection dehned by M and 
a global problem Pq (hL( 7 ); 7 ^), with 1T(7) given by (38), the local problem 
Pl ^{-^*( 7 )}^^ ; 7 ^j, with Xi given by (39), is called a localization of Pq if 
the global admissibility condition in (40) holds. 


Definition 16 (Localization Gap). The localization gap between the global 
problem Pq iW{'y )]and a localization Pl ({-^i( 7 )}^i ; 7 ^) is given by 


\\PL-PG\\=^ll-lh, (46) 

where 7 g and 7 ^ are the respective solutions of the global and local prob¬ 
lems, and are given by 7 g = Pg (hb( 7 ); 7 ^) and jl = Pl ({-^i( 7 )}*^i; 7 ^) = 

max{ 7 j}^p where 7 * = Ph, (Xj( 7 ); 7 ^) is the solution of synthesis problem 
associated with the subsystem Hi. 
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Definition 17 (Localization Distance). Given an interconnection defined by 
M, a global problem Pq ( 1 ^( 7 ); 7 ^) and a localization Pl ^{-^ 1 ( 7 )}^^^ 

the localization distance between the global supply rate W ( 7 ) and the local 
supply rates X = is given by 

D{X, ID) = {QJYl{X)Q, - QJYg{W)Q2) . (47) 

Remark 21 (Closest Localization). The closest localization is then obtained 
by solving the following minimization problem 

min D(X, W) 

s.tQiY^Q.-QjYaQi^O, ^ 

which is equivalent to the semi-definite program 

min —t 
s.t 

tl -QjYLiX)Q^YO (49) 


QJYlQi-QJYgQ2Y0, 

Remark 22 (Exact Localization). The localization is said to be exaet if the 
global admissibility condition holds with equality. In this case the localization 
distance, hence the localization gap, is zero. 

Lemma 1. Let P^(Xb 7 ^) be the closest localization of the global problem 
PciWi^Yn^). Then 

Qj {Yl{X)-Yl{XY)QiYO, (50) 

for all localization Pi(X; 7 ^) of PG(hD( 7 ); 7 ^)- 
Proof. Suppose there exist a localization P/;,(X; 7 ^) such that 

Qj {YL{X)-YL{XY)QiyO. 


It follows that 

Xm^n {QJYl{XYQi) - Xrmn {QJYl{X)Q,) < 0 
which implies that D{X, W) < D{X\ W), a contradiction. 


□ 
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Theorem 2. Given an interconnected system and an associated global prob¬ 
lem. If the closest localization is not feasible, then all feasible local problems 
are not globally admissible. 

Proof. Let be the set of multipliers associated with the closest lo¬ 
calization. For any localization with multipliers it follows from 

Lemma [T] that 




which implies that 

N r 


i=l 


(X,(7)-X;(7)) 


y* 


< 0 , V 7 e 




V(7) 


< 


H T 


A'.' (7) 


, V 7 G M, i = 1,..., iV. 


N 


i=l 


Now, Pl ( \ > ; 7 ^ ) not feasible implies that there exists at least one 


f G {1,..., N} such that 


Yi 

Vi 

9 i 


x!h) 


A.(7) 


Yi 


du < 0, for some Vj G L2e 


du < 0 , 


which implies that Pl 7 -^i( 7 )}^i ; 7 ^j not feasible. □ 

The closest localization has the biggest set of multipliers over which the 
IQCs are dehned. Thus, it’s solution is a global lower bound on the solution of 
all globally admissible local problems. Consequently, the closest localization 
represents the most relaxed local problem that is globally admissible. It 
creates a nice way to transform global control synthesis problem into local 
syntheses problems for the subsystems in such a way that the resulting closed 
loop subsystems will interact through the interconnection to ensure that the 
feasibility of global problem is guaranteed. 
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It is noteworthy that no specihc control design method is enforced. Only 
a somewhat unihed description of of control design problem is made using 
IQC/Dissipativity. This is one of the merits of this approach because it 
allows for heterogeneous control synthesis for the subsystems. Moreover, ad¬ 
dition requirements (e.g robustness, saturation etc) can be built in at the 
subsystems level. The only requirement is the resulting closed loop subsys¬ 
tem satishes the IQC dehned by the corresponding multiplier of the closest 
localization at the respective interconnection input-output port. 

6 .1. Group localization 

Definition 18 (Group). A group is any hnite collection of subsystems 

with group capacity #(G) < N^ where the function ^{.) returns the cardi¬ 
nality of its argument. 

Remark 23. Order does not matter in this dehnition. In other words, two 
groups are considered the same if the elements of any can be produced by 
rearranging the elements of the other. 

Starting with two integers {Ng^N} < N, the objective in this section is 
to identify disjointed groups with 

#(G.) = N. (52, 

and the corresponding multipliers such that the local problem 

Pl ; 7 ^ j is the closest localization. Here, for each localization, 

the jth synthesis problem is dehned for the corresponding group which itself 
is a global problem with respect to its element subsystems. The underlying 
assumption in this section is that the capacity of each group is small enough 
that their respective synthesis problem can be solved without the need for 
any localization. This can be guaranteed by picking N small enough. 

The schematic for the group localization is shown in Figure where 
Pij denotes the degree of membership of subsystem i in group j. For this 
problem, we take pij G {0, l},z G {1, 2,... A^}, j G {1, 2,..., A^}. The 
following constraints are imposed on pij] 



25 



Figure 3: Group Localization 


1. Each subsystem is allowed one and only group assignment. 

= 1, Vze{l,2...,iV}. (53) 

i=i 


2. The number of subsystems in each group cannot exceed the group ca¬ 
pacity N. 


N 

J2P^J<N<N, Vje{l,2...,iV,}. (54) 

i=l 


Proposition 4. Let P = G {0, Pe a symmetric matrix of binary 

values that satishes 


• rank(P) = Ng, 

• P ^ 0, 

• Pa = 1, hi e {1,2, 


( 0 , 


IS a 


Then P admits the factorization P = pp , where p G 
matrix of binary numbers whose elements satisfy the constraints in (53) and 
(54) for some N < N. Moreover, the nonzero singular values of P are 

cr{p) = 
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Remark 24. The elements of matrix P can be interpreted as 


P 




1 if the subsystems i and j belong to the same group 
0 otherwise 


( 55 ) 


Proof. It follows from all the properties of P that there exists a permutation 
matrix T such that 


P 





(56) 


where Xfq. is an Nj x Nj matrix of Is, with Nj = ^(Gj). Equivalently, 


Ia^i 


P = T 






(57) 


where 1 is an iVj-dimensional vector of Is. Thus, the hrst conclusion follows 
by taking 


p = T 






(58) 


Moreover, using (57), it is clear that P is similar to 


’-Ni 




. Thus 


a (P) = a 


/ 

Xni 

\ 

V 

_ 

/ 




(59) 


□ 


The result above demonstrates that the group membership function can 
be inferred from a positive dehnite binary matrix satisfying the states hy¬ 
pothesis. Thus, the group allocation problem is given by 
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min D{PoX, W) 
s.t 

QJYl{P o X)Q^ - QJ Yg{W)Q2 Y 0, 
PYO, 

rank(P) = Ng, 

Pij G {0, 1}, Pii = 1, 

where the binary operator o dehned as 


PoX 


PllXii P^2X^2 ■■■ PinXin 


PinXJj^ P2NXjjg . . . PnnX^n 


(61) 


Notice that the local supply rate is now allowed to be a full block as op¬ 
posed to the diagonal block of the previous sub section, thereby giving more 
room to reduce conservatism. The optimization problem above is nonconvex 
and NP-hard. The details of the relaxation approaches and algorithms to 
solve the problem is left as future work. However, it is noted that, for a 
hxed P, the problem is convex and well-behaved. Thus a quick alternative 
approach considered in this present work is alternating minimization. First, 
the optimization is relaxed using the second part of proposition and then 
regularized as follows; 


min D{PoX, IT) + H(P) 
s.t 

QjY^iP o X)gi - Q^Yg{W)Q 2 y 0, 

P ^ 0, _ ^ ’ 

(^maxiP) < N, 

Pij G [0, 1] , Pjj = 1, 

where f2(.) is any sparsity encouraging norm e.g /i-norm[T7]. trace-norms 
jl8j . etc. Consequently, the alternating minimization algorithm is given in 
Algorithm 




Algorithm 1 Alternating minimization for group localization 
1: Set as such that 

QJYl{X^)Q, - Q^Yg{W)Q2 Y 0 

2: Initialize k = 0 and e > 0 big enough 

3: while e > tol do > tol is a specihed convergence criterion 

4: Set by solving the minimization problem 

min A)(P^oX^hh) + fi(P) 
s.t 

Q7iA(p" o x")gi - qJYg{w)Q2 y 0, 

P^ P 0, 

(TmaxiP^) < N, 

Pt, e [0, 1] , Pt, = 1. 

5: Set by solving the minimization problem 

min P(P^oX^+7hh) 
s.t 

QjYUP’^ o X^+^)Qi - QjYGiW)Q 2 Y 0. 

6 : e^\D{P’^oX^+\W)-D{P^oX’^,W)\ 

7: k ^ k + 1 

8 : end while 
9: return X^ and P^ 
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7. Future Works 


We conclude this paper by giving some interesting directions for future 
works. 

Extended generalization 

The approaches used in this paper can be extended to a more general 
class of problems by replacing the quadratic assumptions to a more general 
convex functions. An interesting example of such generalization is sum of 
squares (SOS). SOS is vastly studied in literature and has very rich applica¬ 
tion in controls [ini EDI ED- As such, there exists a wealth of knowledge to 
build on. Direct application angles, for instance, are; more general descrip¬ 
tion of syntheses problems using integral sum of squares constraint and SOS 
parametrization of the multipliers. 

Components specification for a large-scale design problem 

One can describe components characteristics in a big design project as 
an IQC. Then use the methods developed in this paper to identify the set 
characteristics that best achieve the design objective which is also expressed 
as an IQC. 

Detailed analysis of algorithms 

Detailed algorithm design and analysis for the global admissibility prob¬ 
lem and the group localization problems where not considered. Off-the-shelf 
ADMM and alternating minimization algorithms were used. Design of spe¬ 
cialized algorithms for the problem of the form considered in this paper is an 
interesting problem to consider. 

Application 

Applying a new technique to particular practical problems is always an 
interesting endeavor. The methods presented in this paper describe new ways 
of looking at decentralized control and performance specihcation problems. 
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Specialization to port-controlled Hamiltonian systems 

Port-controlled Hamiltonian Systems (PCHS) are studied extensively [121 
1221 ESI ES] as a generalization of the network modeling of physical systems 
with independent storage elements. Specializing the techniques in this paper 
to PCHS can give very useful insight into the design of control for large 
scale interconnected physical systems. Combined with techniques like graph 
separation [25], the results in this paper can be used to do control synthesis for 
large-scale multi-energy domain physical systems using graphical modeling 
tool like bondgraph[2B] 

Dynamic interconnection 

The problems considered in this paper assumed static interconnection. 
A natural extension is the inclusion of dynamics in the interconnection. It 
is not clear at this point if the results obtain so far will extend to dynamic 
interconnection. Moreover, what happens if the global and local IQCs are 
known and one is interested in obtaining an interconnection (static or dy¬ 
namic) that guarantees the feasibility of the global admissibility condition is 
an interesting question to consider. 

Optimization counterpart of localization 

Given a large scale convex optimization problem. Is there a way to obtain 
a set of smaller optimization problems in such a way that its solution is related 
to the solution of the original problem in a way we can quantify/design? 
If this is true, it might give an interesting extension to operator splitting 
problems. 
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